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Abstract. 

We investigate the dynamics of the extended simple point charge (SPC/E) 
model of water in the supercooled region. The dynamics at negative pres- 
sures show a minimum in the diffusion constant D when the density is 
decreased at constant temperature, complementary to the known maxi- 
mum of D at higher pressures. A similar trend in the rotational diffusion 
is also observed. 



1. Introduction 

The thermodynamic description of supercooled water has been a major 
topic of research already for many years. Most of this scientific effort has 
been concentrated in understanding the static anomalies present in this 
complex fluid. It expands on freezing and, at a pressure of 1 atm, the density 
has a maximum at 4°C. Additionally, there is a minimum of the isothermal 
compressibility at 46° C and a minimum of the isobaric heat capacity at 
35°C [1]. These anomalies are linked with the microscopic structure of liquid 
water, which can be regarded as a transient gel — a highly associated liquid 
with strongly directional hydrogen bonds [2,3]. Each water molecule acts 
as both a donor and an acceptor of bonds, generating a structure that 
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is locally ordered, similar to that of ice, but maintaining the long-range 
disorder typical of liquids. 

Several scenarios have been proposed to explain these anomalies. In 
the " stability limit conjecture " [4,5], the liquid spinodal line for water 
is reentrant. It has a minimum at negative pressures and passes back to 
positive P as the temperature decreases. The increasingly anomalous ther- 
modynamic behavior of liquid water as it is cooled at positive pressures can 
be interpreted in terms of approaching this spinodal. The " critical point 
hypothesis " proposes a second critical point as the terminus of a phase 
coexistence between a high density and a low density liquids. The increase 
in the response functions is then interpreted gnature of being in the 

vicinity of this critical region. The " singularity-free hypothesis " suggests 
that there is no divergence in the response functions. They grow but stay 
finite [3,6,7]. 

Recently, dynamic properties of water have gained attention both exper- 
imentally [8,9] and in computer simulations [10-16]. The surprising result 
of these works is that water also exhibits an anomalous dynamical behavior. 
The increase of the applied pressure leads to an increase in water trans- 
lational diffusion coefficient and to a faster rotational diffusion [8,9,17]. 
These effects can be understood as follows. The increase of pressure leads 
to an increase in the number of defects and in the presence of interstitial 
water. This disrupts the tctrahedral local structure, weakening the hydro- 
gen bonds, and thus increasing the diffusion constant [15,16]. A further 
increase in the pressure leads to steric effects which works in the direction 
of lowering the mobility. The interplay of these factors leads to a maximum 
in the diffusion constant [15, 16] at some high density Pmax- As a result. for 
each isotherm a maximum of the diffusion coefficient and a minimum in the 
rotational correlation time are found. However, the behavior at very low p 
is less well understood. 

In this manuscript, we present our recent studies of how the dynamics 
of low-temperature water are affected by the decrease of the density [18]. 
We perform molecular dynamics (MD) simulations of the SPC/E model of 
water in the range 210 K < r < 280 K and 0.825 g/cm^ < p< 0.95 g/cm^. 
We calculate the rotational and translational diffusion coefficients in this 
region. A relationship between the behavior of the two coefficients is sug- 
gested. 

2. Results 

We performed molecular dynamics simulations using 216 water molecules 
described by the extended simple point charge (SPC/E) model [19], in the 
canonical ensemble (NVT), in a cubic simulation box using periodic bound- 
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ary conditions. The diffusion coefficient D was calculated using the slope of 
the linear regression of the mean square displacement versus time using a 
range of time long enough to assure that the molecules have reached a diffu- 
sional behavior. We show D along isotherms in Figure 1. For T < 260 K, D 
has a minimum value at p 0.9 g/cm^, which becomes more pronounced at 
lower T (Figure 1). This behavior can be understood considering the struc- 
tural changes that occur with decreasing density. At low T, the decreased 
density enhances the local tetrahedral ordering, which leads to a decrease 
in D. Further decreases in density reduces the stability of the tetrahedral 
structure and causes an increase of D. The location of the minimum is near 
the ice Ih density ^ 0.915 g/cm^, which is the density where the perfect 
tetrahedral order occurs. The orientational relaxation was analyzed using 
the rotational autocorrelation functions [20,21]: 

C«(e) = (P,[e(t)-e(0)]) 

where e is a chosen unity vector describing the orientation of the molecules 
and Pi is the i-th order Legendre Polynomial (we restrict ourselves only to 
the study of the first two correlations): 

Pi{x) =x 



^ ' 2 2 

We choose three vectors to describe the orientation. The first vector 
is the unit vector with the same orientation as the dipole moment. The 
second corresponds to the 0-H bond direction and the third is a vector 
perpendicular to the plane of the molecule. 

The correlation functions were fitted to a biexponential decay function [21], 
and two relaxation times were calculated, corresponding to a slower (I) and 
a faster (II) mode. 

C = oo exp(-6t^/2) +a/exp(-t/r^) + an exp{-t/T^^) 

For each correlation using a given vector, two times were calculated, 
leading to a rather complicated symbology. For example, the slower relax- 
ation time for the second-order correlation using the third vector would be 

I (2) 

Tg . With few exceptions, however, we have found only one relaxation 
time, for correlation in a given vector, that is Ty'^^^ = Ty^'^^\ We obtain 
therefore six correlation times for each simulated point, r}^"* and r}^^ for 
the first orientation vector and similarly T2^\ T2^\ and rP. At T = 
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240 K where a detailed analysis was carried out, we find a maximum in the 
orientational time as illustrated in Figure 2. 



3. Conclusions 

Wc analyze the dynamic properties of supercooled water. For high den- 
sities (p > pmax), water behaves as a normal liquid and the decrease of 
translational diffusion coefficient, D, with increasing pressure is governed 
by steric effects. For /Jmin < P < Pmax, as the pressure is decreased, the 
presence of defects and interstitial water decrease, the tetrahedral struc- 
ture dominates, with stronger hydrogen bonds. This process reaches its 
maximum at p = pmin ~ Pice- Further stretching destabilizes the hydrogen 
bond network, leading to an increase in mobility. Preliminary studies of 
the rotational diffusion show a maximum in the orientational time at low 
density region in compass with the behavior of D. 
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5. Caption Figures 

Figure l:(a) Dependence of the diffusion constant D on p along isotherms 
(for p < 1.0 g/cm^). Open symbols arc from ref. [18], and filled symbols 
are from ref. [15]. The dotted line separates liquid state points from phase 
separated state points, but is not an indication of the exact psp{T), which 
varies slightly with T. (b) Full p dependence of D, also showing the maxima. 

Figure 2: Dependence of the rotational relaxations times ri — tq along 
T = 24:0K isotherm as a function of density. 
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